Relaxation time scales in collective dynamics of liquid alkali metals 
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£3 ■ In this paper the investigation of the dynamical processes of liquid alkali metals is executed 

by analyzing the time scales of relaxation processes in liquids. The obtained theoretical dynamic 
£SJ ■ structure factor S(k, u>) for the case of liquid lithium is found to be in excellent agreement with 

the recently received inelastic X-ray scattering data. The comparison and interrelation with other 
theories are given here. Finally, an important part of this paper is the confirmation of the scale 
uniformity of the dynamic processes in liquid alkali metals predicted by some previous molecular 
dynamic simulation studies. 
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I. INTRODUCTION 

The dynamic structure factor S(k,u>) is an experimentally measured term, containing information about the pro- 
cesses in a liquid with long- and short-time scales. It can be used to judge on the microscopic behavior in a system 
\^ . on the basis of its spectra, obtained by means of Inelastic Neutron Scattering (INS)^ or Inelastic X-ray Scattering 
(IXS)£. As for simple liquids, at present a great amount of experimental data of S(k,u>) has been accumulated, in 

• particular, for liquid alkali metals. These data indicate legibly the presence of the collective propagation excitations 
O ' 

beyond the hydrodynamic region. The characteristic feature of liquid alkali metals is a triple-peak shape of S(k,uj) 
lasted to k ~ 0.8k m , where k m corresponds to the first maximum of the static structure factor S(k). Moreover, the 

cannot be obtained within a hydrodynamic treatment, therefore, they are related in some works to the so-called 
kinetic collective excitations. The impossibility to describe these microscopic phenomena and, therefore, to reproduce 
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qualitatively the experimental S(k, uo) by means of ordinary hydrodynamic equations led to the development of other 



theoretical models and approaches. 

One of the simplest and perhaps the earliest modeling approaches is the so-called viscoelastic theory. It allows one to 



obtain the central quasi-elastic line as well as two inelastic peaks symmetrically located around u> = for mesoscopic 
space-frequency region. However, as shown in Refs. 4 '^, this model can not be used for the exact reproduction of 
the experimental spectral shapes of S(k,uj) (see, for instance, the cases of liquid cesium and lithium in Refs4£). 
Therefore, in Ref£ the double-scale model for the viscous relaxation process with fast and slow time scales was tested, 
and as a result a good agreement with the IXS experimental data for the dynamic structure factor was received. 
Recently the similar approach was also applied for the description of relaxation processes in H-bonded liquids®' 7 . The 
existence of two time scales in this model reflects the presence of physically different decay mechanisms. A faster 
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process is hypothetically associated with interactions between an atom and the "cage" of its nearest neighbors, and 
a slower one is identified with the well-known structural (a—) process. However, relaxations of both processes are 
approximated by exponential dependencies. In recent works the viscoelastic model has been improved by means of 
the Markovian closure on the next relaxation level of Zwanzig-Mori hierarchy 8 , it is equivalent to the exponential 
relaxation on this level. It is worth mentioning two others methods, one of which is related to the extension of 
the usual hydrodynamic analytical expressions by modification of hydrodynamic modes to /c-dependence (see, for 
instance, Ref, 9 '). This method assumes the existence of non-hydrodynamical additional modes. The second approach 
is related to the so-called concept of generalized collective modes, which was proposed for the investigation of the 
time correlation functions (TCF's) beyond the hydrodynamic regioniS. The key idea of this method consists in the 
correct choice of the basic set of dynamical variables. 

All these methods are more or less successfully used for the description of collective dynamics in liquids. They have 
common property. Namely, they are actually constructed on heuristic assumption about the presence of exponential 
decay (or combination of exponential decay contributions) in some relaxation processes. Nevertheless, the transition 
and imposition of different relaxation modes in disorder systems can occur even in case of a concrete relaxation 
process, that complicates the selection of the analytical time dependence for the corresponding TCF. This fact is 
proved by the successful application of different mode-coupling theories. On the other hand, this difficulty can be 
resolved by means of analysis and comparison of the resulting time scales of relaxation processes. Therefore, in the 
present work we suggest the approach, which allows us to avoid the immediate approximation of relaxation processes 
by analytical functions. It is based on the development of Bogoliubov's ideas about the hierarchy of relaxation times 
in liquids^, adapted to the formalism of time correlation functions. 

One of the open problems in studying of liquid state (in particular, of the microdynamics of simple liquids) is 
to describe and understand on a general ground the common features of different relaxation processes^. It is well 
known, that the dispersion of the side (high-frequency) peak of dynamic structure factor is the same for all alkali 
metals. Moreover, it is also valid in case of more complex systems, for example, for liquid alloys^. Then the following 
questions arise: Is the origin of relaxation processes the same for liquid systems with the similar features? Can the 
unified description be applied to these systems? As for the group of melting alkalis, it has been indicated in Refi 1 ^ that 
both the equilibrium and the time-dependent correlations can be cast in a properly scaled form for all the the alkali 
metals. Further, it was justified by ab initio molecular dynamic studying in Refi^ too. Experimental confirmation 
of this result was impossible over a long period particularly because of the difficulties related to the technique of 
INS due to the deficient precision of the experimental data. Recently, due to progress in IXS technique this issue 
was considered again 14 . In this work we present investigations related to the determination of corresponding scale 
transitions for liquid systems. 

The organization of the paper is as follows. In the next Section, we describe the theoretical formalism, and the 
comparison with the experimental data and other theories is carried out. The possibility of scale uniformity of 
dynamical processes in the group of liquid alkali metals is analyzed and discussed in Section IIIII The scale-crossing 
relations are also presented here. Finally, we come up with some concluding remarks in Section llVI 
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II. THEORETICAL FORMALISM 



Basic notions 



Let us consider the liquid system of N identical classical particles of the mass m in the volume V and take the 
density fluctuations 

1 N 

WVk) = -=]Te^ (1) 

i=i 

as an initial dynamical variable. To construct a some set of dynamical variables necessary for the description of the 
evolution of the system we use the technique of projection operators of Zwanzig-Morii£*i&. It is a formal version of 
the Gram-Schmidt orthogonalization process, which allows one to obtain the set of orthogonal variables 

W(fc) = {W (k), W 1 (k),W 2 (k), U , .:/,-> } . (2) 

They satisfy the condition (W*Wi) — 5j^i(\Wj\ 2 ) and are connected by the following recurrent relation^: 

W j+1 {k) = CW,{k) - tf^W^k), 

i = 0, 1, 2,...; W- 1 (k) = 0. (3) 

Here the characteristic of the corresponding jth relaxation process, the so-called frequency parameter f2|(fc), appears, 
C is the Liouville operator 

! N V N | 

Tr^r~ E ( 4 ) 
j=l i>j=l ) 

with the momentum of the jth particle and the pair potential u(j, i). So, if Wo(k) is the density fluctuations, then 
Wi(k) is the longitudinal component of the momentum density and so it goes on. 
The TCF's for the corresponding dynamical variables are given by 

(W7(fc)e^V,(fc)) 

M ^ (M) = wimm) 1 Jj= 1 ' 2 '- (5) 

For convenience normalized time correlation functions are used here. The time-evolution operator of Eq. JSJ contains 
the reduced Liouville operator 

£^=(1-^)^(1-^), (6) 



defined by the following projection operators 



Wj(k))(W*(k) 

n ^W' ¥ " w - (7) 

From the condition of orthogonalization of the dynamical variables we obtain the initial values for the TCF's of Eq. 
©: 

>, ifi^, cross-correlations 
M J! (A,* = 0) = ^ . (8) 
1, if j = I, autocorrelations 
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These correlation functions Mji(k,t) are symmetrical in I and j, i.e., 

M j i{k,t) = M lj {k,t). (9) 

Autocorrelation functions of Zwanzig-Mori formalism have the following property: every autocorrelation function of 
the higher order Mj(k,t) = Mjj(k,t) is a memory function for the previous one, i.e. Mj-x(k,t) (autocorrelation 
functions will be marked by one index only in accordance with the used variable), and they are interrelated by 
integro-differential non-Markovian equations of the form: 

dMj ~^ Kt) +^(1,) J* dTM J (k,T)M^ 1 (k,t-T) =0. (10) 

Differentiating the first equation of the chain l|l(J|) . i.e. j — 1, one obtain the generalized Langevin equation: 



^ A y fc ' f ) +nl(k)M (k,t)-nl(k)Ql(k) f dr f dr 'M 2 (k,t-T)M 1 (k,t-T)M (k,T) =0 . (11) 
dt Jo Jq 

One the other hand, these functions describe concrete relaxation processes, the physical meaning of which may be 

established from direct definitions of TCF's. For instance, Mo(k,t) describes the dynamics of fluctuations of density 

correlations in the system, M\(k,t) is the TCF of the fluctuations of the longitudinal component of the momentum 

density, M2(k,t) contains the TCF of fluctuations of energy density. So, these quantities are associated with the 

TCF's of the well-known hydrodynamic "slow" variables. These TCF's have characteristic time scales, which can be 

found from 

/•OO 

Tj(k) = Re / dt Mj(k, t) = Re Mj(k, s = 0), (12) 
Jo 

where Mj(k, s) is the Laplace transform of the corresponding TCF, i.e. Mj(k, s) = J Q °° dt e~ st Mj(k,t)i^i^&. 

So, the memory function approach with single initial dynamical variable extracts the whole set, which describes the 
relaxation processes of the corresponding relaxation levels. In fact, the well-known problem of the choice of a set of 
variables required for the correct description of the system dynamics here is reduced (i) to the search of the number 
of variables for a priori known succession W(fc), that was excellently shown by the recurrent relation approach in a 
works of Le o 21 ' 22 ; and/or (ii) to the finding the correct closure of the chain I|1U|) . 

The ratio between To(k), Ti(fc) and T2(fc) may be quite arbitrary. In the hydrodynamic region (fc — > 0, u — > 0) they 
take large values due to the slow changes of the correspondent variables: densities of mass, momentum and energy. 
Further, one can suggest that the relaxation times of the following TCF's, in comparison with the scales of these 
three variables, are comparable, i.e. r^,(k) ss T^{k). We emphasize here that this assumption does not contradict the 
viscoelastic model, which presupposes that T 2 (k) ^S> ^(k). Obviously, this key condition of the viscoelastic theory is 
just a special case in our approach. Simultaneously, our approach does not deny the presence of the long-lasting time 
tail of M 2 {k, i), which may be adequately taken into account by the mode-coupling theory^. 

Then, taking into account Eq. (|12|) one can find 

M 4 (k,t) = M 3 (k,t) + h(k,t), (13) 

where the "tail" function h(k 7 t) appears. From the short-time asymptotic of the time autocorrelation functions and 
the condition of the long time attenuation of correlation Eq. H13|l yields the following properties of h(k, t): 

lim h(k, t) = lim h(k, t) = 0, (14) 
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this function must have at least one crossing with the time axis at the intermediate region24. Eq. (|13[) allows us to 
obtain the closure of hierarchy of equations of the form 1)1 Oil at the fourth level (j — 4) and by means of Laplace 
transformation to find its exact solution for Mo(k,iuj), in particular, which is directly related to the experimentally 
available term, the dynamic structure factor, S(k,ui). The expression for the resulting S(k,u>) is given in work^i in 
terms of the first four frequency parameters fif(fc), Q 2 (k), f2§(fc), Of(fc) and the Laplace transform of tail function, 
i.e., h(k,iu>). In some cases, the regime with h(k,t) — > may be realized. It can be observed in some parts of time 
(frequency) scale. In this case we find the following expression for the dynamic structure factor: 



s(k,uj) = ^nl(k)n 2 2 (k)nl(k)[4nl(k)-uj 2 ]i{nt(k)ni(k) 

Z7T 

+ uj 2 [ni(k)nl(k) - 2nl(k)nj(k) - nf(k)nt(k) 

+2n 2 1 (k)nl(k)nl{k) - nl(k)nj(k)nl(k) + n%(k)nl(k)] 
+ uj 4 [nl(k) - 2nl(k)nl(k) + 20?(fc)^(fc) - 2n 2 2 (k)nl(k) + nl(k)nl(k)] 

+ u> 6 [fll(k) - SlKk)]}- 1 . (15) 

This equation is also expressed through the first four frequency parameters, which are directly related to the first five 
even frequency moments of dynamics structure factor. It is necessary to note that this expression is obtained in the 
way completely different from the theory of moments^. 



B. Comparison with IXS experiment and relationship with other theoretical approaches 

In Fig. 1 we report the dynamic structure factor S(k,uj) of liquid lithium (T = 475K) for some wave numbers 
calculated from Eq. l|15fl (solid line) and obtained from IXS experiment (circles)—. Being used in theoretical com- 
putations the static structure factor S(k) for both cases was taken from Ref^£. The first frequency parameter was 
directly defined from its definition f^(fc) = KsTk 2 /mS{k). The second frequency parameter 0|(fc) is related to the 
fourth frequency moment. We found this parameter from the values of the infinite frequency sound velocity Coo(fc)^*^ 
by means of relation Coo(fc) = yfVl\(k) + Sl\{k)/k. The high-order parameters were found by comparison with the 
experiment. Eventually, we have revealed that all the frequency parameters have the similar dispersion. In particular, 
they have the first principal maximum at the same wave numbers such as the side peak of S(k, lu), i.e. at k ~ 0.55/c m , 
and any low order parameter is less than the high order one. 

We would like to emphasize that the theoretical S(k,ui) and, in particular, the position of the side peak, is very 
sensitive to the magnitude of fi|(fc). The magnitudes of Q 2 {k) and fl 2 (k) influence the form of S(k,u}). However, it 
is not so important to know these parameters separately as their ratio, i.e. ri|(fc)/il|(fc). 

To compare the theoretical outcome with the experiment we modified it to account for the quantum mechanical 
detailed balance condition according to 

and then broadened it for the finite experimental resolution effects R(k, uj)^: 

J R(k,uj- u')S q (k,Lj')du. (17) 
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From Fig. 1 one can see that the above described theoretical approach yields a good agreement with IXS data of 
both systems. 

Now we can execute a more detailed study of the obtained results and compare them with other approaches: the 
usual viscoelastic model, the double-scale viscous model and the generalized mode approach. The common feature of 
these theories is the use of the time autocorrelation function Mji(k, t) of Eq. at j = I = 2. So, the viscoelastic and 
the double- viscosity models are based on approximations to this term, and M 2 (k, t) plays a key role in these theories. 
As for our approach, it gives the following form for Laplace transform of M 2 (k, t) 



M 2 {k,s) = [s + CllityM^s)]- 1 
s + nl{k)M 3 {k,s) 



(18a) 



which are obtained by Laplace transform of the third and fourth (j — 3, 4) equations of the chain (|10() . 

To pass from the frequency dependence of M 2 (k,iuj) to the time one, let us consider the low frequency region 
restricted by the value 2D,^{k). For convenience we introduce here a small parameter (at the fixed wave number k): 

e = ^, d9) 



Taking into account the fact that the found values of f2|(fc) for liquid sodium and lithium achieve 10 29 — 10 30 s~ 2 for 
the low-fc region, we span by introducing parameter £ the frequency (time) range uj < 10 15 s _1 (t > 10~ 15 s), which is 
important for us and is available experimentally. 

Expanding the radicand in Eq. I|18b|) as a series in the parameter £ 



v / T+e = i + |-y + ■■■ , (20) 



we can rewrite it in the following way 



M3(s) = -2^ + i + M-^ + -' (21) 
By restricting the number of terms in the series (|2U[) [and, accordingly, in Eq. |JHJ] we receive from Eq. (|18a|l the 
linear combination of the Lorentz functions 

. A (k) 

M 2 (k,s) = J2 Zl I77T> 3 = 1,2,3,5,..., (22) 

j S + T J ( fc ) 

the number of which will be increased at the increase of the number of terms in the series (120(1 . The quantities Aj(k) 
and Tj(k) are expressed by the relaxation frequencies ^ 2 (fc) and Vl\(k). Going over to the time scale by the inverse 
Laplace transform 2 ^ we obtain 

M 2 (k,t) = J2 A j( k ) e ~ t/TAk) - (23) 

3 

By restricting the first term of the series l|2(J|l only we receive the simplest model from the first equality of Eq. I|18a|) 
with Eq. 

M 2 (k,t) = e-*/ T « (24) 
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which corresponds to the viscoelastic model with the relaxation time r(fc) = Q,4,(k)/D^(k), and from the second equality 
of Eq. Ijl8a(l the double exponential model, i.e Eq. (|23JI at j = 2, with the following time relaxation parameters 



n,2(A) = 

and the weight factor 



(25) 



n A (k) + ^n'i(k) -n 2 3 (k) 



A(k) = " w v 4W J W . (26) 

(a) - ni(fc) 

This case may be related to the double-time viscous model 5 *^, two-time exponential ansatd^^. In the general 
form Eq. I|23[) corresponds to the framework of generalized collective mode approaches with the sum of the weighed 
exponents for the TCF M%(k, t), where r7 (fc) denote eigenvalues of a generalized dynamic matrix with the elements 
consisting of static correlation functions, and the weight factors Aj(k) are the amplitudes describing the contribution 
of the corresponding modes. 

So, it is obvious that the theory underlying Eq. Ijl5|) prescribes such behavior of the second order memory function 
M%{k, t), which may be represented in the form of Eq. (|23|) and can be reduced to the above-mentioned models. Eq. 
(|23JI is in fact an expansion of M2(k,t) into decay channels embedded in this function. 



III. SCALE UNIFORMITY OF DYNAMICS PROCESSES IN LIQUID ALKALI METALS 

The determination of the scale uniformity of structural and dynamical features for different groups of liquids is very 
important for the physics of liquid state. On the one hand, it allows one to apply the unified theoretical description 
to the whole group. On the other hand, it allows one to remove the difficulties related to obtaining the experimental 
data. The fact is that until recently the microscopic dynamics of liquids could be experimentally probed by INS only. 
However, there were often different problems related to, firstly, separation of collective and one-particle contributions, 
and, secondly, gross experimental errors (and even with impossibility to obtain data) for different (fc, to) -regions. 
Recent progress in the technique of IXS has allowed one to clear some of the obstacles^. Ten years ago the possibility 
of the unified description of the structural and dynamical properties of different liquid alkali metals near the melting 
point was found by the comprehensive molecular dynamics simulation study^, where the adopted potential model of 
Price, Singwi and Tosi was used, and the scale passage was executed on the basis of the potential parameters. The 
recent sketchy attempt of testing this outcome experimentally has shown its inconsistencjii*. 

In present work we also execute the comparison of the dynamic structure factor spectra of liquid lithium and sodium. 
As known from the experimental results, the dynamic structure factor S(k,u>) depends strongly on the temperature 
T and the wave number k. So, one can define the reduced forms of these terms as T/T m and k/k m , where T m is the 
melting temperature and k m is the main peak position in the static structure factor S(k) for the corresponding system. 
The scale time interval t* can be expressed as t* — fc _1 \Jm,j KbT. Thus defined time unit t* is different from the 
one introduced in Ref. 14 , because the present term varies with the change of space and temperature characteristics. 
Although we do not exclude the possibility, that this scale unit may be independent of the temperature and the wave 
number for other systems (for instance, semi-conductors, or H-bonded liquids). 

In Fig. 2 we report the comparison of S(k,ui) spectra for liquid lithium and sodium 514 at approximately the same 
reduced temperatures T/T m and wave numbers k/k m . Namely, T/T m — 1.049 for liquid lithium and 1.051 in case of 
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sodium. From this figure one can see that dynamic structure factor practically coincides in the first two higher cases. 
From the lower plot of Fig. 2 one can see, that the position of inelastic and central peaks for both systems is the same. 
However, though the overall coincidence of spectra is observed at intermediate frequencies only, the peak altitudes 
are a little different. Such deviation can easily be explained by the fact that the plot for liquid lithium is presented 
for a higher value of the reduced wave number, 0.75, whereas in case of sodium k/k m — 0.73. As known, these wave 
numbers correspond to the so-called the dc Gennes narrowing region characterized by a strong A:— dependence. In 
other words, a higher section at k of flat S(k,u))/t* is presented for lithium than for sodium. It is necessary to take 
into account, that the values of the reduced temperatures for both systems are also slightly different. 

Notice that the time unit t* depends on the system features (m), on the probed spatial region (k) and the tempera- 
ture regime (T) in contrast to scale units k m and T m , which remain unchanged the spatial region and the temperature 
of the system is revised. 

At result, the experimental or theoretical S(k,u>) for any single metal allows one easily to restore this term for the 
whole group of alkali metals at same reduced conditions, k/k m and T/T m . Moreover, the theory developed for the 
concrete separate alkali metal may be simply extended to the whole group. 

As an example, in Fig. 3 we report the dynamic structure factor of liquid potassium S K (k,iv) obtained on the 
basis IXS data for liquid sodium S Na (k, oj). The transition S Na (k, u) — ► S K (k, ui) has been executed by means of the 
following scale reductions: 



h.Na \ rn K r rNa 

S K (k,.) = S^(k,u) jr] J- Fni¥F , (27a) 



K _ , Na K I 171 1 
UJ — CO 



f,Na y m K T Na ' 



(27b) 



T K = ** = (27c) 

By the top subscript we note the corresponding system (K or Na). 

IV. CONCLUDING REMARKS 

The following results are presented in this work. 

(i) The theory, developed on the basis of Bogoliubov ideas about the hierarchy of relaxation times, allows one 
to obtain dynamic structure factor, reproducing adequately experimental IXS spectra for liquid alkali metals (in 
particular, for liquid lithium and sodium) in the region of low values of wave number. 

(ii) The expansion of the second order memory function into exponential decay channels, used (sometimes intu- 
itively) in others theories, may be easily obtained within the framework of the presented approach. This is the evidence 
of the multi-mode character of decay of the observed relaxation process. 

(iii) An important result of this work is the confirmation of the proposition about the unitary description of the 
dynamical features of liquid alkali metals, and finding of corresponding scale transition relations. 
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VI. FIGURE CAPTIONS 

Fig. 1. Dynamic structure factor of liquid lithium at the temperature T = 475K. The solid lines are the results of 
the theoretic model (|15|l , whereas the open circles are the IXS data^. The theoretical lineshapes have been modified to 
account for the quantum mechanical detailed balance condition and broadened for the finite experimental resolution 
effects as described in the text. The wave numbers k are given in a reduced form, where k m is the main peak position 
in the static structure factor S(k). 

Fig. 2. IXS spectra of liquid lithium at T = 475K (□□□) and liquid sodium at T = 390K (o o o)Ui in the reduced 
units. The scale frequency to* is chosen as the term inverse proportional to t* . 

Fig. 3. The dynamic structure factor S(k, us) of liquid potassium at T = 354. IK calculated from IXS data of liquid 
sodium at T = 390K by the scale reduction described in the text. 
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